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Letter to the Editor: 
Convective envelopes in rotating OB stars 
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ABSTRACT 

Aims. We study the effects of rotation on the outer convective zones of massive stars. 

Methods. We examine the effects of rotation on the thermal gradient and on the Solberg-Hoiland term by analytical 
developments and by numerical models. 

Results. Writing the criterion for convection in rotating envelopes, we show that the effects of rotation on the thermal 
gradient are much larger and of opposite sign to the effect of the Solberg-Hoiland criterion. On the whole, rotation 
favors convection in stellar envelopes at the equator and to a smaller extent at the poles. In a rotating 20 Mq star at 
94% of the critical angular velocity, there are two convective envelopes, with the bigger one having a thickness of 13.2% 
of the equatorial radius. In the non-rotating model, the corresponding convective zone has a thickness of only 4.6% of 
the radius. The occurrence of outer convection in massive stars has many consequences. 
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1. Introduction 

It is generally considered that the Cowling model applies 
to massive OB stars: i.e., a convective core surrounded by 
a large radiative envelope. However, long since stellar mod- 
els have shown that massive stars have an outer convec- 
tive en velope encom passing several perce nt of the stellar 
radius (|MaederJll980l ). Also. iLangerl (|l997l ) has shown that 
an Eddington factor V — kL/(AttcGM) tending toward 1.0 
implies convection. Our aim is to show that fast rotation 
amplifies the size of the convective envelope in OB stars as 
well as to develop anisotropic convective envelopes. 

Various limits can be considered about the ef fects of 
rotat i on and high luminosity on the stellar stability (jLangerl 
Il997t iMaeder fc MevnetJ l2000t) : the T-limit, which is the 
Eddington limit for T — > 1; the O-limit, which is reached by 
stars at rotational break-up with a small or negligible effect 
of the Eddington factor T; the SIT- limit, which applies to 
stars where both luminosity and rotation play significant 
roles. We show that not only the stars at the T— limit, but 
also the stars at the Or-limit and at the O-limit, have 
amplified external convective zones. 

The occurrence of outer convective envelopes in OB 
stars and their anisotropic structure lead to many astro- 
physical consequences: 

— Convection generates acoustic modes that may allow 
asteroseismic observations of OB stars. 

— Convective motions may play a role in driving mass loss 
by stellar winds. 

— For stars close to the critical rotation, convective mo- 
tions lower the effective break-up velocities. 



— An outer convective envelope may make a dynamo and 
contribute to some chromospheric activity generating 
an X-ray emission from OB-stars. 

— A convective envelope transports chemical elements and 
angular momentum. 

— The occurrence of outer convectio n may modify the von 
Zeipel theorem (|von Zeipell[T92l . 

A closer investigation is justified. We start by an analytical 
approach (Sect. [2]), and finish by two-dimensional models 
of 20 M Q rotating envelopes (Sect. [3]). 



2. Convection in rotating stars 

In a rotating star of mass M, luminosity L, and angular 
velocity SI (supposed to be shellular, i.e. constant on shells), 
the total gravity is the sum of the gravitational, centrifugal, 
and radiative accelerations: 



9tot — 9ef[ + 9rad — Sgrav + Srot + ffrad • 



(1) 



The vector g e fi has both radial and tangential components, 
the radial component at colatitude $ is 



9c«, 



GM r 



Q 2 r 3 2 



(2) 



where r is the radius at colatitude The radiative accel- 
eration is directed outward 
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where F is the flux. On an isobaric surfa ce, F is given by 
the von Zeipel theorem (jvon Zeipell 119241 ) . 



L(P) 



with M*(r) = M r 1 - 



4ttGM* 

n 2 



9eff , 



2nGp I1 



(4) 
(5) 



and L(P) is the luminosity on the isobar, p m the inter- 
nal av erage density . In baroclinic stars, there are other 
terms (jMacdcr 1990), however they are small and neglected 
here. The flux is proportional to the effective gravity g e ff- 
The effective mass M* is the mass reduced by the cen- 
trifugal force. Let us note that one has VL 2 / '(2nGp m ) « 
(4/9)(w/w crit ) 2 ~ (16/81)w 2 , where w = il/fl crit , v is 
the rotation velocity at the level considered and i> cr it = 
(2/3)[GM/i* crit (tf = 0)]V2. 

2.1. Effect of rotation on the thermal gradient 

Formally the Solberg-Hoiland criterion is to be considered 
in a rotating star, as in Sect. 12.21 However, the radiative 
gradient V ra d is also modified by rotation, an effect gen- 
erally not accounted for in Schwarzschild's criterion. The 
local flux and the equation of hydrostatic equilibrium are 



F=- X VT and VP=gg cS , 



(G) 



with x — 4acT 3 / (3kq). Radiation pressure is included in 
P, the total pressure. The local radiative gradient becomes 
in a rotating star, 



_dT dnP _ 
Vrad = ~An~~dPT = 



3 kL(P)P 
16nacG M*(r)T 4 ' 



(7) 



where the derivatives are computed along a direction n per- 
pendicular to the isobars. Except L(P), the terms are local 
and thus have to be taken at a given (r, ■d). We ignore the 
horizontal thermal gradient and take L(P) as constant in 
the envelope. With ([5]) and the expression of the Eddington 
factor T, we get 



V ra d 



4(1-/3) (l 



o 2 

2irGp„ 



(8) 



where (3 = P g /P is the ratio of gas pressure to the total 
pressure, thus P/(aT 4 ) = 1/[3(1 — /?)]. The adiabatic gra- 
dient V a d is 



'ad 



6/3 



32 - 24/3 - 3/3 2 



(9) 



As T varies with (3 also varies with colatitude, and we 
write That /?(#) is higher at the equator favors equa- 

torial convection. The criterion for convective instability 
V r ad > Vad becomes 



i - 



o 2 



> 4[l-/3(i?)]V 



ad 



(10) 



where the ^-dependence of T comes only through k($) 
(|Maeder fc MevnetH2000fh Eq. (TO) has various interesting 
consequences: 
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Fig. 1. The ratio (3 = P gas /P as a function of the stellar 
masses on the zero-age main-sequence with Z — 0.02. The 
continuous line shows the value at the stellar centers, the 
dotted line the maximum value of /3 inside these models. 
The long-dashed line indicates the minimum value of (3 
permitted by equilibrium conditions (|Chandrasekhar|[l984t 
!MitalaslfT997h . 



— In the absence of r otation, expre ssion © is equivalent 
to Langer's result (|Langerlll997l ). The right-hand side 
of Eq. pop is always smaller than 1.0, thus if T — > 1, 
the criterion is satisfied. Convection is present in layers 
close to the Eddington limit. 

— In a rotating star, inequality (|10|) is more easily satisfied. 
Thus, rotation favors convection in stellar envelopes. 

— The occurrence of convection depends on both k($) 
and f3(i}). Equatorial ejection is always favored, even 
for electron scattering opacity caused by the higher (3. 

— When the centrifugal force can be derived from a po- 
tential (conservative case) , the temperature and density 
are constant on isobars and so that T and (3 are also 
constant on isobars. In that case, rotation also favors 
convection as can be seen from Eq. (9). 

One can wonder whether rotating stars that are neither 
at the critical nor at the Eddington limit may develop a 
convective envelope. The ratio (3 decreases with mass (Fig. 
[IJ, while the Eddington factor T increases, e.g. T = 2.5 x 
10" 5 , 0.0047, 0.021, 0.098, 0.239, 0.343, 0.544 for 1, 5, 9, 20, 
40, 60, and 120 M Q stars on the ZAMS. The parameter (3 
is at a minimum in the stellar centers, reaches a maximum 
in the envelope, and is zero at the stellar surface. One has 
the following relation between the maximum /3-values in 
the outer layers and T, 







1.0-sT 



1 



(1-/3) 



(11) 



with s = 0.72 ± 0.01 between 20 M and 120 M . Using 
Eq. (fTTj) . one eliminates T from criterion (flO|) and gets 



1 



1 



o 2 

2irGp„ 



> 4 



6/3 



32 - 24/3 - 3/3 2 



(12) 



This relation indicates above which value of u> there is con- 
vection for a given value of (3 at the maximum (for lower 
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(3 in the envelope, the inequality is evidently satisfied more 
easily). For example, for (3 — * 1.0, the inequality becomes 

2 ^Gp > 0-132 or (w/wcrit) > 0.54. This is an approxima- 
tion owing to the simplified relation adopted and because 
the various parameters also vary with depth. However, it 
shows that rotating massive stars not even at the critical 
limit may have enhanced convection. 



2.2. The Solberg-Hoiland criterion 

A fluid element displaced in a rotating star is also sub- 
ject to the restoring effect of angular momentum conserva- 
tion . This leads to the Solberg-Ho iland criterion for stabil- 
ity (|Kippenhahn fc Weigertl [l990T ). which is (for constant 
mean molecular weight fi) 

Vad - Vrad + Vfj Sin 1?>0 (13) 



with Vr 



Hp i d(n 2 w 4 ) 



(14) 



5grav<5 dw 

where w — r sin -d is the distance to the rotation axis and 
5 = — (9 In q/ dh\T) p . The quantity Vn depends on the 
distribution of the specific angular momentum j = nj 2 Sl, 
which results from transport processes. As j decreases out- 
ward, Vn generally has a stabilizing effect. Let us consider 
the two extreme cases for Sl(r): 

- 1. Constant specific angular momentum: a distribution 



SI 



may result from the Rayleigh- Taylor instability. 



This distribution is also sometimes considered in convective 
regions, with the argument that the plumes rapidly redis- 
tribute the angular momentum. If so, Vn = 0, and one is 
brought back to Schwarzschild's criterion. 
-2. Constant angular velocity: this assumption is also used 
in convective regions, with the argument that turbulent vis- 
cosity favors solid rotation. If so, Vn simplifies to 

„ A n 2 Hp 4 si 2 p 

Vn = 4 f = . (15) 

ffgrav ffgrav Q 5effO 

We can simplify this expression further. In the outer layers, 
as long as k w const, and g c g w const, at an optical depth 
r one has P ss (ffeff/ '«) t. This gives 

si 2 fn 2 R 3 " 

Vn 



GM 



Q K RS 



(16) 



Sgrav 9 K5 

The term in the first parenthesis is to 2 , while that in square 
brackets is just the ratio (R — r)/R (assuming S = 1), which 
is small in the envelope. The Solberg-Hoiland criterion be- 
comes in this approximation, 

R 



1 



Q 2 



27rGp„ 



> 4 (1 - 0) Vad + 



R 



sin , (17) 



where as above the various quantities are local ones. At low 
rotation, the Solberg-Hoiland term Vn is negligible with 
respect to the other terms. At high rotation for constant 
SI, it is not negligible, but in general smaller than the other 
terms because the convective zone lies very close to the 
surface and the term (R — r/R) is small. 

Since the actual rotation laws are likely between the 
two extreme cases f2(r) = const, and SI ~ r -2 , we conclude 
that the main effect of rotation on convection in stellar 
envelopes is not the inhibiting effect due to the Solberg- 
Hoiland criterion, but the effect of rotation on the thermal 
gradient (Eq. [TT)|l. which enhances convection. 



3. Numerical models 

We do some 2-D models of the outer regions of a 20 M Q 
fast-rotating star with X = 0.70 and Z = 0.020 (Fig. HJ. 
At each latitude we integrate the equations of the structure 
for the corresponding effective gravity and T c s of the Roche 
model of the given rotation, also taking the effect of the 
reduced mass into account. In the envelope, we suppose 
that is a constant as a function of depth (the problem is 
conservative). We first consider only the effect of rotation 
on the thermal gradient and then the complete Solberg- 
Hoiland criterion to see the differences. 



3.1. Effects of rotation on the thermal gradient 

Without rotation, a 20 Mq model at the end of the MS 
evolution has two outer convective zones. The first one is 
very close to the surface and is due to an increase of the 
opacity caused by partial He ionization. It extends from 
r/R = 0.992 to 0.999, i.e. only 0.7% of the radius, and 
contains a very small fraction of the total stellar mass (2.5 ■ 
10 -9 ). The second one is deeper, between r/R = 0.915 to 
0.962 (4.7% of the total radius), and contains a fraction 7.4- 
10~ 7 of the total mass. Both convective zones are associated 
with opacity enhancements. 

For fast rotation, with a ratio Sl/S7 cr it = 0.94, where 
Sl cr it is the critical angular velocity, the convective lay- 
ers are shown in Fig. [21 They are more extended than 
without rotation. The thin upper convective zone extends 
from r/R = 0.987 to 0.999 at the pole, i.e., over 1.2% of 
the stellar radius, and it contains 1.3 • 10~ 8 of the stel- 
lar mass. The deeper convective zone covers the region be- 
tween r/R = 0.836 to 0.936 (i.e. 10.0% of the polar ra- 
dius), and its mass fraction is 2.8 • 10~ 6 . At the equator, 
we have the following sizes for the two convective layers: 
between r/R = 0.958 and r/R = 0.988 for the first one 
(3.0% of the equatorial radius) and between r/R = 0.727 
and r/R = 0.862 for the deeper one (13.5% of the radius). 
The included masses are the same as at the pole. 

In agreement with Sect. 12. 1[ convection is more ex- 
tended in the rotating model than in the non-rotating one. 
Fig-E]shows that contrarily to the classical Cowling model, 
massive rotating stars have a large-size convective enve- 
lope. Interesting is that, if we look at the structure of the 
envelope of the star as a function of the pressure, we see that 
this structure is independent of the colatitude. This is ex- 
pected since we have supposed SI constant in the envelope. 
In that recalled above, T and (3 are constant on iso- 

bars and the extensions of the convective zones, expressed 
in term of differences of pressure between the bottom and 
the top of the convective zone are the same at the pole 
and at the equator. The spatial extensions are, however, 
greater in the equatorial region than in polar ones. This 
comes from the variations in the spatial gradients of the 
pressure and temperature with the colatitude imposed by 
the hydrostatic equilibrium (gradients of pressure have to 
balance pg c s)- Another result of the constancy of pressure 
and temperature (and thus density) on an isobar is that the 
radiative gradient is also constant on this isobar, except for 
the change in the effective mass as given by M», which is 
lower than M in rotating stars. Thus, for rotating stars the 
radiative gradient is larger and convection is favored not 
only at the equator, but at each colatitude, compared with 
the non-rotating model. 
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Fig. 2. 2-D model of the external convective zones and of 
the convective core (dark areas) in a model of 20 M Q with 
X = 0.70 and Z = 0.020 at the end of MS evolution with 
fast rotation ( oj = f2/f2 crrt = 0.94). The axes are in units 
of cm. 

The mass loss rate in the considered model is 6.2 x 
10~ 7 Mq yr" 1 . This means that within one year about 4 
times all the matter in the thin convective zone is lost in 
the stellar winds! Thus, the matter carried by the winds 
is continuously passing through the superficial convective 
zones in a dynamical process. 

3.2. Solberg-Hoiland criterion 

We also computed the envelope structure with the Solberg- 
Hoiland criterion for convection. The values of (V a d — V ra( j) 
and Vnsini? are shown in Fig. [3] for a 20 M Q model at 
the equator, i.e. where the effect of the Solberg-Hoiland 
criterion is the strongest. The solid line shows the dif- 
ference between the adiabatic and the radiative gradient 
(identical to the values obtained in the model computed 
with the Schwarzschild criterion). The short-dashed line 
shows the Solberg-Hoiland term. The limits of the con- 
vective zones when the Soilberg-Hoiland criterion is used 
are inside the regions where (V a d — V ra d) is negative. 
Indeed, the limits are where (V a d — V ra d) + Vn situ? = 0. 
Convective zones are thus smaller than those obtained with 
the Schwarzschild criterion. We find the following values for 
the extension of the convective zones at the equator when 
the Solberg-Hoiland criterion is used: between r/R = 0.960 
and r/R — 0.988 for the superficial one (2.8% of the equa- 
torial radius) and between r/R = 0.727 and r/R = 0.859 
for the deeper one (13.2% of the total radius), i.e. values 
slightly smaller than but very similar to those obtained in 
the model computed with the Schwarzschild criterion. This 
numerical example confirms that the Solberg-Hoiland term 
Vn sin i? has a very limited influence in stellar envelopes, as 
discussed in Sect. 2.3. 

4. Conclusions 
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Fig. 3. V a d — Vrad (solid line) and Vnsini? (short-dashed 
line) at the equator in a 20 M Q model with f2/Q cr j t = 0.94 
at the end of the MS. The long-dashed line indicates the 
zero level. 

sign with respect to the Solberg-Hoiland criterion, so that 
rotation favors convection instead of inhibiting it. The in- 
crease of the convective zone occurs mainly at the equator 
and also a bit at the poles. In a fast-rotating 20 M Q Pop I 
star, there are two equatorial zones covering a total of 16% 
of the stellar radius at the equator. 

There are several consequences of thees results to be ex- 
amined in future. The outer convective motions may lower 
the escape velocity as well as the critical rotation veloc- 
ity. The matter accelerated in the winds continuously goes 
through the convective zone in a dynamical process, sug- 
gesting that convection plays a role in accelerating the stel- 
lar winds and in producting the clumps in the winds. The 
convective pistons generate acoustic waves of periods of sev- 
eral hours to a few days. The density is very low, and it is 
thus likely that convection injects oscillations into the wind 
rather than into the interior. 
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In stellar envelope of rotating stars, the effects of rotation 
on the thermal gradient arestronger and with the opposite 



